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5_! ' Abstract 

In this paper, we search for absolute limitations of the Tile Assembly Model (TAM), along 
with techniques to work around such limitations. Specifically, we investigate the self-assembly 



of fractal shapes in the TAM. We prove that no self-similar fractal fully weakly self-assembles 
at temperature 1, and that certain kinds of self-similar fractals do not strictly self- assemble at 
any temperature. Additionally, we extend the fiber construction from Lathrop et. al. (2007) to 
show that any self-similar fractal belonging to a particular class of "nice" self-similar fractals 
' has a fibered version that strictly self-assembles in the TAM. 

C0 ! 



1 Introduction 



Self-assembly is a bottom-up process by which (usually a small number of) fundamental components 
(NJ | automatically coalesce to form a target structure. In 1998, Winfree [18] introduced the (abstract) 

Tile Assembly Model (TAM) - an extension of Wang tiling \16\ I17j . and a mathematical model 
of the DNA self-assembly pioneered by Seeman et. al. [2]. In the TAM, the fundamental com- 
^ | ponents are un-rotatable, but translatable "tile types" whose sides are labeled with glue "colors" 

and "strengths." Two tiles that are placed next to each other interact if the glue colors on their 
abutting sides match, and they bind if the strength on their abutting sides matches, and is at least 
a certain "temperature." Rothemund and Winfree (13^ 1 12j later refined the model, and Lathrop et. 
al. [10] gave a treatment of the TAM in which equal status is bestowed upon the self-assembly of 
^ . infinite and finite structures. There are also several generalizations [21 [TT| [8] of the TAM. 

Despite its deliberate over-simplification, the TAM is a computationally and geometrically ex- 



pressive model. For instance, Winfree [18] proved that the TAM is computationally universal, and 
thus can be directed algorithmically. Winfree [18] also exhibited a seven-tile-type self-assembly 
system, directed by a clever XOR-like algorithm, that "paints" a picture of a well-known shape, 
the discrete Sierpinski triangle S, onto the first quadrant. Note that the underlying shapes of each 
of the previous results are actually infinite canvases that completely cover the first quadrant, onto 
which computationally interesting shapes are painted (i.e., full weak self-assembly). Moreover, 
Lathrop et. al |9] recently gave a new characterization of the computably enumerable sets in terms 
of weak self-assembly using a "ray construction." It is natural to ask the question: How expressive 
is the TAM with respect to the self-assembly of a particular, possibly infinite shape, and nothing 
else (i.e., strict self-assembly)? 
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In the case of strict self-assembly of finite shapes, the TAM certainly remains an interesting 
model, so long as the size (tile complexity) of the assembly system is required to be "small" relative 
to the shape that it ultimately produces. For instance, Rothemund and Winfree [13] proved that 
there are small tile sets in which large squares self-assemble. Moreover, Soloveichik and Winfree [15] 
established the remarkable fact that, if one is not concerned with the scale of an "algorithmically 
describable" finite shape, then there is always a small tile set in which the shape self-assembles. 
Note that if the tile complexity of an assembly system is unbounded, then every finite shape trivially 
(but perhaps not feasibly) self- assembles. 

When the tile complexity of an assembly system is unbounded (yet finite), only infinite shapes 
are of interest. In the case of strict self-assembly of infinite shapes, the power of the TAM has only 
recently been investigated. Lathrop et. al. [10] established that self-similar tree shapes do not 
strictly self-assemble in the TAM given any finite number of tile types. A "fiber construction" is 
also given in [10], which strictly self-assembles a non-trivial fractal structure. 

In this paper, we search for (1) absolute limitations of the TAM, with respect to the strict 
self-assembly of shapes, and (2) techniques that allow one to "work around" such limitations. 
Specifically, we investigate the strict self-assembly of fractal shapes in the TAM. We prove three 
main results: two negative and one positive. Our first negative (i.e., impossibility) result says that 
no self-similar fractal fully weakly self-assembles in the TAM (at temperature 1). In our second 
impossibility result, we exhibit a class of discrete self-similar fractals, to which the standard discrete 
Sierpinski triangle belongs, that do not strictly self-assemble in the TAM (at any temperature). 
Finally, in our positive result, we use simple modified counters to extend the fiber construction 
from Lathrop et. al. [TO] to a particular class of discrete self-similar fractals. 

2 Preliminaries 

2.1 Notation and Terminology 

We work in the discrete Euclidean plane 1? = ZxZ. We write Ui for the set of all unit vectors, i.e., 
vectors of length 1, in Z 2 . We regard the four elements of U2 as (names of the cardinal) directions 
in I?. 

We write [X] 2 for the set of all 2-element subsets of a set X. All graphs here are undirected 
graphs, i.e., ordered pairs G = (V,E), where V is the set of vertices and E C [V] 2 is the set of 
edges. A cut of a graph G = (V, E) is a partition C = {Cq,C\) of V into two nonempty, disjoint 
subsets Co and C\. 

A binding function on a graph G = (V, E) is a function (3 : E — ► N. (Intuitively, if {u, v} £ E, 
then (3 ({it, v}) is the strength with which u is bound to v by {u, v} according to 0. If is a binding 
function on a graph G = (V, E) and C = (Co, C\) is a cut of G, then the binding strength of on 
C is 

c = {0(e)\eeE,enC o ^^ andenCi/0}. 
The binding strength of (3 on the graph G is then 

f3{G) = min{/? c |C is a cut of G} . 

A binding graph is an ordered triple G = (V,E,0), where (V,E) is a graph and (3 is a binding 
function on (V, E). If r £ N, then a binding graph G = (V, E, (3) is T-stable if 0(V, E) > r. 



A grid graph is a graph G = (V, E) in which V Q Z 2 and every edge {rh, n} G E has the 
property that rh — h G The /uZZ grid graph on a set C Z 2 is the graph Gy = (V, -E) in which 
E contains every {rh, n} G [V] 2 such that m — n G L^- 

We say that / is a partial function from a set X to a set Y , and we write / : X —-> Y, if 
/ : D — > Y for some set DCI. In this case, D is the domain of /, and we write D = dom /. 

All logarithms here are base-2. 



2.2 The Tile Assembly Model 

We review the basic ideas of the Tile Assembly Model. Our development largely follows that of 
|12| . but some of our terminology and notation are specifically tailored to our objectives. In 
particular, our version of the model only uses nonnegative "glue strengths", and it bestows equal 
status on finite and infinite assemblies. We emphasize that the results in this section have been 
known for years, e.g., they appear, with proofs, in 



Definition. A tile type over an alphabet £ is a function t : U2 — > £* X N. We write t = (cob, str^), 
where cob : U2 — > S*, and str t : U2 — > N are defined by t(u) = (colt(w), strait)) for all u £ Ui- 

Intuitively, a tile of type t is a unit square. It can be translated but not rotated, so it has a 
well-defined "side u " for each u € Ui- Each side u of the tile is covered with a "glue" of color 
cob (u) and strength str^ (u) . If tiles of types t and t' are placed with their centers at rh and rh + u, 
respectively, where rh € Z 2 and u S U2, then they will bind with strength str^(u) • {t(u) = t'(—u)J 
where [</>] is the Boolean value of the statement <f>. Note that this binding strength is unless the 
adjoining sides have glues of both the same color and the same strength. 

For the remainder of this section, unless otherwise specified, T is an arbitrary set of tile types, 
and r G N is the "temperature." 

Definition. A T- configuration is a partial function a : Z 2 — » T. 

Intuitively, a configuration is an assignment a in which a tile of type a(m) has been placed 
(with its center) at each point rh G dom a. The following data structure characterizes how these 
tiles are bound to one another. 

Definition. The binding graph of a T-configuration a : Z 2 — » T is the binding graph G a = 
(V,E,{3), where (V,E) is the grid graph given by 

V = dom a, 

E = I {m, n} G [y] 2 m — n G C/ n , col a (^) (n — rh) = col Q (^) (rh — h) , and str Q (^) (n — rh) > o| , 

and the binding function (3 : E — > Z + is given by 

/3 ({m, n}) = str Q{ ^) (n - m) 

for all {rh, n} G E. 
Definition. 



1. A T-configuration a is r- stable if its binding graph G Q is r-stable. 



2. A r-T- assembly is a T-configuration that is r-stable. We write A? for the set of all r-T- 
assemblies. 

Definition. Let a and a' be T-configurations. 

1. a is a subconfiguration of a', and we write a C a', if dom a C dom a' and, for all m G dom a, 
a(m) = a!(m). 

2. a' is a single-tile extension of a if a C a' and dom a' — dom a is a singleton set. In this case, 
we write a' = a + (rh \— > t), where {rh} = dom a' — dom a and £ = a'{m). 

Note that the expression a + (rh i— ► i) is only defined when rh El? — dom a. 
We next define the "r-i-frontier" of a r-T-assembly a to be the set of all positions at which a 
tile of type t can be "r-stably added" to the assembly a. 

Definition. Let a G .Ay. 

1. For each i G T, the T-t-frontier of a is the set 



3[a 



m G Z 2 — dom a 
2. The t- frontier of a is the set 



strt(u) • [a(m + u)(—u) = t(u)j > r 



9 T a = (9^Q. 



The following lemma shows that the definition of d\a achieves the desired effect. 

Lemma 1. Let a G A T , m G Z 2 — dom a, and i G T. Then a + (m i— > i) G A T if and only if 
rh G 9^"a. 

Notation. We write a — > a' (or, when r and T are clear from context, a — !-> a') to indicate 

T,T 

that a, a' G A T and a' is a single-tile extension of a. 

In general, self-assembly occurs with tiles adsorbing nondeterministically and asynchronously 
to a growing assembly. We now define assembly sequences, which are particular "execution traces" 
of how this might occur. 

Definition. A r-T-assembly sequence is a sequence a = (a^ | < i < k) in A r , where k G Z + U{oo} 
and, for each i with 1 < i + 1 < k, an — — » aj + i. 

Note that assembly sequences may be finite or infinite in length. Note also that, in any r-T- 
assembly sequence a = (ati \ < i < k), we have a, C aj for all < % < j < k. 

Definition. The result of a r-T-assembly sequence a = (q^ | < i < k) is the unique T- 
configuration a = res(a) satisfying dom a = Uo<j<fc dom a.{ and ajCa for each < i < k. 

It is clear that res(a) G A T for every r-T-assembly sequence a. 



Definition. Let 

1. A t-T- assembly sequence from a to a' is a r-T-assembly sequence a = (oij | < i < k) such 
that ao = a and res(a) = a'. 

2. We write a ► a' (or, when r and T are clear from context, a — ► a') to indicate that there 

exists a r-T-assembly sequence from a to a'. 

Definition. An assembly a € A T is terminal if and only if d T a = 0. 

We now define tile assembly systems. 

Definition. 

1. A generalized tile assembly system (GTAS) is an ordered triple 

T = (T,a,r), 

where T is a set of tile types, a £ is the seed assembly, and r G N is the temperature. 

2. A taZe assembly system (TAS) is a GTAS T = (T, <r, r) in which the sets T and dom cr are 
finite. 

Intuitively, a "run" of a GTAS T = (T, a, r) is any r-T-assembly sequence a = (a>i \ < i < k) 
that begins with ao = a. Accordingly, we define the following sets. 

Definition. Let T = (T, a, r) be a GTAS. 

1. The set of assemblies produced by T is 

A[T] = [ae A? 



a > a 

2. The set of terminal assemblies produced by T is 

= {a G -4,[T]| a is terminal} . 

Definition. A GTAS T = (T, cr, r) is directed if \Aa[T\\ = 1. 

We are using the terminology of the mathematical theory of relations here. The reader is 
cautioned that the term "directed" has also been used for a different, more specialized notion in 
self-assembly [?]. 

In the present paper, we are primarily interested in the self-assembly of sets. 
Definition. Let T = (T, a, r) be a GTAS, and let X C I?. 

1. The set X weakly self-assembles in T if there is a set B C T such that, for all a € »4o[T], 
a^OB) =X. 

2. The set X fully weakly self-assembles in T if X and Z 2 — X both weakly self-assemble. 



3. The set X strictly self-assembles in T if, for all a E .*4rj[T], dom a = X. 



Intuitively, a set X weakly self-assembles in T if there is a designated set B of "black" tile types 
such that every terminal assembly of T "paints the set X - and only the set X - black" . Moreover, 
a set fully weakly self-assembles in T if every terminal assembly of T tiles the entire plane on top 
of which X is "painted." In contrast, a set X strictly self-assembles in T if every terminal assembly 
of T has tiles on the set X and only on the set X. Clearly, every set that strictly self- assembles in 
a GTAS T also weakly self-assembles in T. 

We now have the machinery to say what it means for a set in the discrete Euclidean plane to 
self-assemble in either the fully weak, weak or the strict sense. 

Definition. Let ICZ 2 . 

1. The set X weakly self-assembles if there is a TAS T such that X weakly self-assembles in T. 

2. The set X fully weakly self- assembles if there is a TAS T such that X fully weakly self- 
assembles in T. 

3. The set X strictly self-assembles if there is a TAS T such that X strictly self-assembles in T. 
Note that T is required to be a TAS, i.e., finite, in all three parts of the above definition. 

2.3 Local Determinism 

The proof of our construction uses the local determinism method of Soloveichik and Winfree |15j . 
which we now review. 

Notation. For each T-configuration a, each rh E Z 2 , and each u E U2, 

str Q (m, u) = stT a r^\(u) ■ \a{fh){u) = a(rh + u)(— u)J. 
(The Boolean value on the right is if {m, rh + u] ^ dom a.) 

Notation. If a = (a^ | < i < k) is a r-T-assembly sequence and rh € Z 2 , then the a-index of rh 
is 

i$(rh) = min{i E N | rh E dom a« }. 
Observation 1. rh E dom res(a) < 00. 

Notation. If a = (ati \ < i < k) is a r-T-assembly sequence, then, for fh,rh' E Z 2 , 

rh rh' <3> is(m) < i^(rh'). 

Definition. (Soloveichik and Winfree [15]) Let a = (c*i \ < i < k) be a r-T-assembly sequence, 
and let a = res (a). For each location rh E dom a, define the following sets of directions. 



1. IN a (m) = j-u E U2 rh + u rh and str^. (m, u) > 1. 

2. OUT"(m) = /« E J7 2 -«6lN s (m + «)}. 



Intuitively, IN a (m) is the set of sides on which the tile at m initially binds in the assembly 
sequence a, and OUT" (to) is the set of sides on which this tile propagates information to future 
tiles. 

Note that IN" (to) = for all rh G a . 

Notation. If a = (oj | < i < k) is a r-T-assembly sequence, a = res(a), and to, G dom a — 
dom ao, then 

a\rh = a \ ^dom a — {to} — (^rh + OUT a (TO)^ . 

(Note that a \ m is a T-configuration that may or may not be a r-T-assembly 

Definition. (Soloveichik and Winfree [E]). A r-T-assembly sequence a = («j | < i < k) with 
result a is locally deterministic if it has the following three properties. 

1. For all m, G dom a — dom ao, 

ueIN s (m) 

2. For all to G dom a — dom ao and all t G T — {a(m,)}, m ^ c^(<3 \ 

3. d T a = 0. 

That is, a is locally deterministic if (1) each tile added in a "just barely" binds to the assembly; 
(2) if a tile of type to a t a location rh and its immediate "OUT-neighbors" are deleted from the 
result of a, then no tile of type t ^ to can attach itself to the thus-obtained configuration at location 
to; and (3) the result of a is terminal. 

Definition. A GTAS T = (T, a, r) is locally deterministic if there exists a locally determinstic 
r-T-assembly sequence a = (Qj | < i < k) with ao = o~. 

Theorem 1. (Soloveichik and Winfree [15]) Every locally deterministic GTAS is directed. 

2.4 Discrete Self-Similar Fractals 

In this subsection we introduce discrete self-similar fractals. 

Definition. Let 1 < c G N, and ICff 2 (we do not consider N 2 to be a self-similar fractal). We 
say that X is a c-discrete self-similar fractal, if there is a set \ i G {0, . . . ,c — 1}} ^ V C 

{0, . . . , c - 1} x {0, . . . , c - 1} such that 

oo 

x = [jx t , 

i=0 

where Xi is the i th stage satisfying Ao = {(0,0)}, and Aj + i = X{ U (Aj + c l V). In this case, we 
say that V generates X. A is a discrete self-similar fractal if it is a c-discrete self-similar fractal 
for some c G N. 

In this paper, we are concerned with the following class of self-similar fractals. 

Definition. A nice discrete self- similar fractal is a discrete self-similar fractal such that ({0, . . . , c— 
1} x {0}) U ({0} x {0, . . . , c - 1}) C V, and G* is connected. 



□ 
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(a) Nice 



(b) Non- 



Figure 1: The first stages of discrete self-similar fractals. The fractals in (a) are nice, whereas (b) shows 
two non-nice fractals. 



2.5 Zeta-Dimension 

The most commonly used dimension for discrete fractals is zeta-dimension, which we use in this 
paper. The discrete-continuous correspondence mentioned in the introduction preserves dimension 
somewhat generally. Thus, for example, the zeta-dimension of the discrete Sierpinski triangle is 
the same as the Hausdorff dimension of the continuous Sierpinski triangle. 

Zeta-dimension has been re-discovered several times by researchers in various helds over the past 
few decades, but its origins actually lie in Euler's (real-valued predecessor of the Riemann) zeta- 
function [6] and Dirichlet series. For each set ACZ 2 , define the A- zeta- function Qa '■ [0, oo) — > [0, oo] 



by Ca(s) = £ 



(0,0)^(m,n)e4.' 



\m\ + \n\) s for all s G [0, oo). Then the zeta-dimension of A is 
Dim^(A) = inf{s [ (a(s) < oo}. 



It is clear that < Dim^(A) < 2 for all iCZ 2 , It is also easy to see (and was proven by Cahen 
in 1894; see also [H[7j) that zeta-dimension admits the "entropy characterization" 



Dim^(A) 



lim sup 

n— >oo 



log \A 



<n| 



logn 



(2.1) 



where A< n = {(k, I) G A \ \k\ + \l\ < n}. Various properties of zeta-dimension, along with extensive 
historical citations, appear in the recent paper [5], but our technical arguments here can be followed 
without reference to this material. We use the fact, verifiable by routine calculation, that (2.1) can 
be transformed by changes of variable up to exponential, e.g., 



Dim c (A) 



lim sup 

n— >oo 



log |A [0 . 



2 n lnNl 



also holds. 



3 Impossibility Results 

In this section, we explore the theoretical limitations of the Tile Assembly Model with respect to 
the self-assembly of fractal shapes. First, we establish that no discrete self-similar fractal fully 
weakly self-assembles at temperature r = 1. Second, we exhibit a class C of discrete self-similar 
fractals, and prove that if F G C, then F does not strictly self-assemble in the TAM. 

Definition. (Lathrop et. al. [10]) Let G = (V,E) be a graph, and let D C V. For each r G V, the 
D-r-rooted subgraph of G is the graph Go,r = (^D,r 5 Ez> r ), where 

Vd,t = {v EV \ every simple path from v to (any vertex in) D in G goes through r } 

and£ D ,r = En[V D>r ] ■ B is a D-subgraph of G if it is a -D-r-rooted subgraph of G for some r G V. 



Definition. Let G = (V, E) be a graph. Fix a set D C V, and let r, r' € V. 

1. (Adleman et. al. PQ) G£> jr is isomorphic to Guy, and we write G\d i7 . ~ G£). r ' if there exists 
a vector a € Z 2 such that V^.r = r' + 

2. We say that Gi) jr is unique if, for all r' € V, Gd )T . ~ G^y =4> r = r'. 

We will use the following technical result to prove that no self-similar fractal weakly self- 
assembles at temperature r = 1. 

Lemma 2. (Adleman et. al. p]) Let X C N such that G\ is a finite tree, and assume that 
X strictly self-assembles in the TAS T = (T,o~,t). Let a € «4rj[7~]. If a (it) = a(v), then the 

The following construction says that if it is possible to self-assemble a finite path P at tempera- 
ture 1 (not necessarily uniquely), then there is always a TAS Tp in which P uniquely self-assembles 
at temperature 1. 

Construction 1. Let T be a finite set of tile types, and a = (pci \ < i < k) be a 1-T-assembly 
sequence, with a = res(a), satisfying 

1. dom ao = {(0,0)}, and 

2. Ga is a connected, finite path P. 

It is clear that for all v £ P, lN a (v) = OXJT a (v) = 1. Now define, for each v € P, the (unique) 

vectors Vm,v ut, satisfying v- m G IN Q ('u), and v ou t £ OXJT a (v). For each v E P, define the tile type 
ttf, where for all u 6 f/2, 

{(col a (#) (u) ■ "in" , str a( fl) (n)) if w + n = v in 

(col Q (y) (n) • "out" , str a( y) (u)) if v + n = v out 

(A, 0) otherwise. 

Let Tp = {ty I u G P}. Note that since P is finite, so too is Tp. Now define the TAS Tp = 
(Tp,ap, 1), where for all v £ N 2 , <7p is defined as 

[I otherwise. 

It is routine to verify that Tp is directed (i.e., P uniquely self-assembles in Tp). 

We now have the machinery to prove our first impossibility result. 

Theorem 2. If F C N 2 is a discrete self-similar fractal, is connected, and F fully weakly 
self-assembles in the TAS Tp = (T, cr, r), where cr consists of a single tile placed at the origin, then 
r > 1. 

Proof. Suppose that P is generated by the set V C {0, . . . c — l} 2 , and assume for the sake of 
obtaining a contradiction that r = 1. Let V' = {0, . . . c — l} 2 — V. There are two cases to consider. 



Case 1 If there exists a path P = ((xq, yo), . . . , (a?z_i, Vi-i)) in Gy,, with Gp connected, satisfying 
either of the following. 

1- Oo, yo) e ({0} x {0, ... c - 1}) and {xi- X ,Ul-x) G ({c - 1} x {0, . . . c - 1}). 
2. (x , yo) G ({0 . . . , c — 1} x {0}) and (a:z_i,ift_i) G ({0, . . . , c - 1} x {c - 1}). 

Without loss of generality, assume that P satisfies (1). First note that there exists a £ V, 
and there is no path from (0,0) to a in Gy. Define, for all i G N, the points 

Si = c l ■ a. 

Since F is infinite, it is possible to choose k G N large enough so that the path 

P = ((xq, yo), (xi,yi), ■ ■ ■ (xk-i,yk-i)) 
satisfies the following properties. 

1. (x ,y ) = (0,0), 

2. there exists I £ N such that (xfc-i, y&-i) = 

3. Gp is connected and simple (in fact a tree), and 

4. there exists a sub-path P' C P, such that G*, is connected, P' C N 2 -P, and |P'| > 12|T| 
(because F fully weakly self-assembles). 

Since r = 1, there is an assembly sequence a = (a, | < i < k), with a = res(a), satisfying 
ao = <t, and dom a = P. Then by Construction Q] there exists a 1-Tp-assembly sequence 
«P = (an | < i < k), with result ap = res(ap) satisfying dom ap = P, and ap(x;_i, yn) G 
P. By (4), there exist s,t G P' such that ap(s) = ap(i), and s, f G" P. Let Pd om o-,s and 
P domo .^be dom u-subgraphs of P. Then Lemma [2] tells us that Pjomo-,? ~ ^dom o-t*' whence 
there exists a location 6 G P' such that ap(6) G P. This contradicts the definition of P. 

Case 2 If there is no such path in Gy,, then we proceed as follows. First note that there exists 
a . It is clear that, for all i G N, c l ■ a + (1, 1) F. For each i G N, define the point 

a i = c l -a + (1,1). 

Since F is infinite, it is possible to choose k G N large enough so that the path 

P = ({x , y ), (xi, y-i), . . . (x k -i,yk-i)) 
satisfies the following properties. 

1. (x Q ,yo) = (0,0), 

2. there exists I G N such that (xfe_i, yu-i) = &1 (because F fully weakly self-assembles), 

3. Gp is connected and simple (in fact a tree), and 

4. for all u G U 2 , min{i | i ■ u + a t G P} > 12|T|. 



Since r = 1, there is an assembly sequence a = (cti | < i < k), with a = res(a), satisfying 
Qo = cr, and dom a = P. Then by Construction [1] there exists a 1-Tp-assembly sequence 
5p = (aj | < i < fc), with result ap = res(dp) satisfying dom ap = P. By (4), there 
exist s,t G P such that ap(s) = ap(i), and s, f G" F. Let Pdom a,s an d -P dom CT f be dom cr- 
subgraphs of P. Then Lemma [2] tells us that P can be extended to an infinite, periodic path 
P' consisting of all but finitely many non-black tiles (i.e., tiles that are placed on the points 
in N 2 — F). This contradicts the definition of F. 

□ 

Note that Theorem [3] says that even if one is allowed to place a tile at every location in the first 
quadrant, it is still impossible for self-similar fractals to weakly self-assemble at temperature 1. 

Next, we exhibit a class C of (non-tree) "pinch-point" discrete self-similar fractals that do not 
strictly self-assemble. Before we do so, we establish the following lower bound. 

Lemma 3. If X C 1? strictly self-assembles in the TAS T = (T, a, r), where a consists of a single 
tile placed at the origin, then |T| > 



| B B is a unique dom cr-subgraph of G* X 



Proof. Assume the hypothesis, and let a G A\j[T]. For the purpose of obtaining a contradiction, 



suppose that \T\ < 



. By the Pigeonhole Principle, 



|s B is a unique dom cr-subgraph of X 

there exists points r, f € X satisfying (1) a[f) = a(r 1 '), and (2) G^ om a p ^ G^ oma ^. Let a' be 
the assembly with dom a 1 = {r 1 }, and for all u G U2, define 



(col a (^/)(u),str Q ( ? ,)(u)) if f + u G G dom a ^ 
nrJW \ (A,0) otherwise. 

Let a" be the assembly with dom a" = {r"}, and for all u G U2, define 



- / { col a(r")(u),str a ^,^(u)) if r" + u G G dom a j„ 
° [T >[U) X ( A '°) otherwise. 

Then T = (T, cr, r) is a TAS in which Gd om o-,r' strictly self-assembles, and T" = (T, a" ,t) is a TAS 
in which Gdom a,f" strictly self-assembles. But this is impossible because a (r 1 ) = a {r") implies 
that, for all u EU 2 , a' (/) (u) = a" (r") (u). □ 

Our lower bound is not as tight as possible, but it applies to a general class of fractals. Our 
second impossibility result is the following. 

Theorem 3. If X C N 2 is a discrete self-similar fractal satisfying (1) {(0, 0), (0, c— 1), (c — 1, 0)} C 
V, (2) V n ({1, . . . c - 1} x {c - 1}) = 0, (3) V n ({c - 1} x {1, . . . , c - 1}) = 0, and (4) G* is 
connected, then X does not strictly self-assemble in the Tile Assembly Model. 

Proof. By Lemma [3l it suffices to show that, for any m G N, 



D 



B is a unique dom cr-subgraph of G* j 



> m. 



Define the points, for all k G N, = c (c(c — 1), c — 1), and let 



B k = {(a,b)£F 



(a,b) G{0,...c fc -l} 2 +f fc } 



Conditions (1), (2), and (3) tell us that Gf is a dom u-subgraph of G* (rooted at r&), and it is 
routine to verify that, for all k, k! £ N such that A; 7^ k', G* 9^ Gg ^ . Thus, we have 



< k < m 



< 



Ib B is a unique dom u-subgraph of G* | 



□ 

Corollary 1 (Lathrop, et. al. [ID]). The standard discrete Sierpinski triangle S does not strictly 
self-assemble in the Tile Assembly Model. 



4 Every Nice Self-Similar Fractal Has a Fibered Version 

In this section, given a nice c-discrete self-similar fractal X C N 2 (generated by V), we define its 
fibered counterpart X. Intuitively, X is nearly identical to X, but each successive stage of X is 
slightly thicker than the equivalent stage of X (see Figure [2] for an example). Our objective is to 
define sets Fq,F\,. . . C Z 2 , sets Tq,T\, . . . C Z 2 , and functions I, f, t : N — ► N with the following 
meanings. 

1. Tj is the i th stage of our construction of the fibered version of X. 

2. Fi is the fiber associated with Tj. It is the smallest set whose union with Tj has a vertical left 
edge and a horizontal bottom edge, together with one additional layer added to these two 
now-straight edges. 

3. is the length (number of tiles in) the left (or bottom) edge of Tj U Fi. 
4- /(») = \F\. 

5. t(i) = \Ti\. 

These five entities are defined recursively by the equations 
To = X2 (the third stage of X), 

F = ({-1} x {-1, . . . , c 2 }) U ({-1, . . . , c 2 } x {-!}) , 
1(0) = c 2 + 1, /(0) = 2c 2 + 1, i(0) = (\V\ + l) 2 , 
T i+1 =TiU ((Ti U Fi) + l(i)V) , 

Fi+i = Fi U ({-» - 2} x {— f — 2, —i — 1, - - ■ ,Z(t + l)-i-3}) 

U ({-» - 2, -i - 1, • • • , i(i + 1) - i - 3} x {-i - 2}) , 

l(i + l)=c-l(i) + l, 

f(i + l)=f(i) + c-l(i + l)-l, 

t(t + i) = |y|t(») + /(»)• 
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Figure 2: Construction of the fibered Sierpinski carpet. The blue, and orange tiles represent (possibly 
translated copies of) Fq, and F\, respectively. Note that this image should be viewed in color. 



Finally, we let 



x = Qt. 



i=0 



Note that the set Tj U Tj is the union of an "outer framework," with an "internal structure." 
One can view the outer framework of Tj U Tj as the union of a square <5j (of size % + 2) , a rectangle Xi 
(of height i + 2 and width — (i + 2)), and a rectangle (of width i + 2 and height Z(i) — (i + 2)). 
Moreover, one can show that the internal structure of Tj U Tj is simply the union of (appropriately- 
translated copies) of smaller and smaller Xj and lj-rectangles. 

We have the following "similarity" between X and X. 

Lemma 4. If X C N 2 is a nice self-similar fractal, then Dim^(X) = Dim^(X). 

In the next section we sketch a proof that the fibered version of every nice self-similar fractal 
strictly self-assembles. 

5 Sketch of Main Construction 

Our second main theorem says that the fibered version of every nice self-similar fractal strictly 
self-assembles in the Tile Assembly Model (regardless of whether the latter strictly self-assembles). 

Theorem 4. For every nice self-similar fractal IcN 2 , there is a directed TAS in which X strictly 
self-assembles. 

We now give a brief sketch of our construction of the singly-seeded TAS Tx = (Xx,«r, 2) in 
which X strictly self-assembles. The full construction is implemented in C++, and is available at 



the following URL: http://www.cs.iastate.edu/~lnsa 



Throughout our discussion, S#, Yg, and X^ refer to the square, the vertical rectangle and the 
horizontal rectangle, respectively, that form the "outer framework" of the set ((Tj U Tj) + ■ u) 
(See the right-most image in Figure H]). 



5.1 Construction Phase 1 

Here, directed graphs are considered. Let X be a nice (c-discrete) self-similar fractal generated by 
V . We first compute a directed spanning tree B = (V, E) of Gy using a breadth- first search, and 
then compute the graph B R = (V, F R ) , where 

E R = {(v,u) | (u, v) E E and u 7^ (0, 0)} U {((0, 1), (0, c — 1)), ((1, 0), (c — 1), 0)}. 
Figure [3] depicts phase 1 of our construction for a particular nice self-similar fractal. 
Notation. For all 7^ u £ V, Ui n is the unique location v satisfying (u, v) £ F R . 
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Figure 3: Phase 1 of our construction. Notice the two special cases (right-most image) in which we define 
(0,l)i„ and (l,0) in . 



5.2 Construction Phase 2 

In the second phase we construct, for each (0, 0) 7^ u £ V, a finite set of tile types that 
self-assemble a particular subset of X. There are two cases to consider. 

Case 1 In the first case, we generate, for each u G V — {(0, 0), (0, 1), (1, 0)}, three sets of tile 
types Ts s , Tx a , and Ty a that, when combined together, and assuming the presence of 
((Tj U Fi) + ■ -Uin), self-assemble the set ((Tj U F{) + ■ u), for any i EN. 

Case 2 In the second case, we generate, for each u £ {(0, 1), (1, 0)}, the same three sets of tile 
types (Ts a , Tx a , and Ty a ) that self-assemble the set ((Tj U Fi) + ■ u) "on top of" the set 
((Ti_i U + i(« ~ 1) ■ ^in), for any » 6 N. 

Finally, we let T x = U(o,o)^«Gy T «> wnere r « = ^ U Tx a U T Ya - Figure H gives a visual 
interpretation of the second phase of our construction. Our TAS is 7x = (Tx,«r, 2), where a 
consists of a single "seed" tile type placed at the origin. Our full construction yields a tile set of 
5983 tile types for the fractal generated by the points in the left-most image in Figure [H 

5.3 Details of Construction 

Note that in our construction, the self-assembly of the sub-structures S^, Y$, and can proceed 
either forward (away from the axes) or backward (toward the axes). 
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Figure 4: Let V be the left-most image. The first arrow represents phase 2 of the construction. The 
second arrow shows a magnified view of a particular point in V. Each point (0,0) ^ u £ V can be viewed 
conceptually as three components: the tile sets Ts fl , Tx u and Ty a that ultimately self-assemble the square 
S^, and the horizontal and vertical rectangles and Y$ respectively. 



times, where p(j) is the largest number of consecutive least-significant 
0's in the base-c representation of j, and [</>] is the Boolean value of the 
statement (j). The value of a row is the number that it represents. We 
refer to any row whose value is a multiple of c as a spacing row. All 
other rows are count rows. The type of the counter that self-assembles 
is u. 

Each counter self-assembles on top (or to the right) of the square S^, 
with the width of the counter being determined by that of the square. It 
is easy to verify that if the width of is i + 2, then TV? self- assembles 
a rectangle having a width of i + 2 and a height of 



5.3.1 Forward Growth 

We now discuss the self-assembly of the set ((Tj U F{) + u ■ l(i)) for u G V satisfying 
(t?+{(-l,0),(0,-l)}). 

If u & {(0,0), (0, 1), (1,0)} (i.e., case 1 of phase 2), then the tile set 
Ts a self-assembles the square fig directly on top (or to the right) of, 
and having the same width (height) as, the rectangle Yq. (X^. n ). If 
u G {(0, 1), (1, 0)} (i.e., case 2 of phase 2), then the tile set self- 
assembles the square on top (or to the right) of the set Y% in such that 
right (top) edge of the former is flush with that of the latter. Note that 
in case 2, the width of Y^ in is always one less than that of S^. In either 
case, it is straightforward to construct such a tile set Ts a . 

The tile set of Ty a self-assembles a fixed-width base-c counter (based 
on the "optimal" binary counter presented in [5j) that, assuming a width 
of i G N, implements the following counting scheme: Count each positive 
integer j, satisfying 1 < j < & — 1, in order but count each number 
exactly 

[c divides j] • p(j) + [c does not divide jj ■ 1 



(c 2 + 1) c 4 + 



1 



l(i)-(t + 2), 



Figure 5: Example of a 
base-3 modified binary 
counter. The darker 
shaded rows are the 
spacing rows. 



which is exactly Y^. Figure 5 shows the counting scheme of a base-3 counter of width 3. We 
construct the set Tx a by simply reflecting the tile types in Ty a about the line y = x, whence 
the three sets of tile types Tg_, Tx g , and Ty a self-assemble the "outer framework" of the set 
((Ti UFi) + u-l(i)). 

The "internal structure" of the set ((Tj U Fi) + u ■ l(i)) self-assembles as follows. Oppositely 
oriented counters attach to the right side of each contiguous group of spacing rows in the counter 
(of type it) that self-assembles Y$. The number of such spacing rows determines the height of 
the horizontal counter, and its type is (0, j/c mod c), where j is the value of the spacing rows to 
which it attaches. We also hard code the glues along the right side of each non-spacing row to 
self-assemble the internal structure of the points in the set To. 

The situation for is similar (i.e., a reflection of its vertical counterpart), with the exception 
that the glues along the top of each non-spacing row are configured differently than they were for 
Yfi. This is because nice self-similar fractals need not be symmetric. 

One can prove that, by recursively attaching smaller oppositely-oriented counters (of the ap- 
propriate type) to larger counters in the above manner, the internal structure of ((Tj UFj) + u- l(i)) 
self-assembles. 



5.3.2 Reverse Growth 

We now discuss the self-assembly of the set ((Tj U Fi) + u ■ l(i)), for all u 6 V satisfying u- m € 
(u + {(1,0), (0,1)}). 

In this case, the tile set Ty a (Tx a ) self-assembles the set Yq (Xjj) directly below (or to the left 
of) the square Sq. , and grows toward the x-axis (or y-axis) according to the base-c counting scheme 
outlined above. We also configure Ty tt (Tx a ) so that the right (or top)-most edge of [Xtf) is 
essentially the "mirror" image of its forward growing counterpart (See Figure [6]). This last step 
ensures that the internal structure of ((Tj U Fi) + u- l(i)) self-assembles correctly. Next, the square 
Sn attaches to the bottom (or left)-most edge of Y$ (X^). Finally, the set X^ (Y$) self-assembles 
via forward growth from the left (or top) of the square S^. 
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(a) 



(b) 



(c) 



Figure 6: (a) depicts forward growth, (b) shows what happens if the tile set Tx s were to simply "count in 
reverse," and (c) is the desired result. 



5.3.3 Proof of Correctness 

To prove the correctness of our construction, we use a local determinism argument. The details of 
the proof are technical, and therefore omitted from this version of the paper. 



6 Conclusion 



In this paper, we (1) established two new absolute limitations of the TAM, and (2) showed that 
fibered versions of "nice" self-similar fractals strictly self-assemble. Our impossibility results moti- 
vate the following question: Is there a discrete self-similar fractal ICN 2 that strictly self-assembles 
in the TAM? Moreover, our positive result leads us to ask: If X C N 2 is a discrete self-similar frac- 
tal, then is it always the case that X has a "fibered" version X that strictly self-assembles, and 
that is similar to X in some reasonable sense? 
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